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Abstract. A search for the rare decays B s — > y + y~ and B — > p. + pT performed in pp collisions at = 
7 TeV is presented. The data sample, collected by the CMS experiment at the LHC, corresponds to an integrated 
luminosity of 1.14 ftr 1 . In both cases the number of events observed after all selection requirements is consistent 
with expectations from background and standard model signal predictions. The resulting upper limits on the 
branching fractions are BF(B S -> fi + iT) < 1.9 x 1(T 8 and BF(B -> n+yT) < 4.6 x 1(T 9 at 95% confidence level 
(CL). Furthermore, the combination of the results of the search for the decay B s — > p + p~ by the CMS and LHCb 
experiments is presented. The combined upper limit is BF < 1.1 x 10~ 8 at 95% CL. 



1 Introduction 

In the standard model (SM) of particle physics, flavor chang- 
ing neutral current decays are highly suppressed (Fig. [TJ: 
they are forbidden at tree level and can only proceed through 
higher-order loop diagrams, they are helicity suppressed 
by factors of (ra//m#) 2 , where m/ and mg are the masses 
of the lepton and B meson and they require an internal 
quark annihilation within the B meson. The SM predictions 
(Tab.[T} are significantly enhanced in several extensions of 
the SM, although in some cases the decay rates are low- 
ered. A (blind) search for the rare decays B s — > yt yT and 
Bq — > ytyC is presented here, using 1.14 ftr 1 of integrated 
luminosity collected by the CMS experiment in 201 1 fl~). 



2 A blind analysis 

We perform a counting experiment in the dimuon mass dis- 
tribution centered on the B s {Bq) meson mass (see Tab.[2]for 
the region definitions). The background is estimated from 
the sidebands and from Monte Carlo (MC) simulation. 
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Fig. 1. Illustration of the rare decays — > y + y~. In the SM 
these decays proceed through box (left) and penguin (right) inter- 
actions. 

Table 1. Branching fraction SM predictions for B s — > y + yT and 
So -> ti+fT. 

Decay channel BF SM predictions 

B s ^>y + y~ (3.2 ± 0.2) x KT 9 
B ^p + p- (1.1 ± 0.1) X 10" 10 



Table 2. Region definitions 
Region definition Invariant mass (GeV) 

overall window 4.90 < m mm < 5.90 

blinding window 5.20 <n lim < 5.45 

B s — » p + pT window 5.20 < m lim < 5.30 

B a — > y + y~ window 5.30 < m^ m < 5.45 



The analysis uses a relative normalization to the well- 
measured decay B ± — > J/iffQj.+fTjK* to avoid a depen- 
dence on the uncertainties of the bb production cross sec- 
tion and luminosity measurements. The B s — > n + n" branch- 
ing fraction is measured using 

N<! f e B± 

BF(B S ^ jtV) = -J r J JL2"BF(B t ^ J/i/tK*) (1) 

N ots f* *»>< 

where BF(B ± — > J/ifr^) is the branching fraction B ± — > 
J/tf/^; Ns/N i, s is the background-subtracted number of 
observed signal candidates in the signal window over the 
number of reconstructed B ± events; ef ot je tot is the ratio of 
the signal total efficiency over the B ± efficiency and f u /f s 
is the ratio of the B ± and B s meson production fractions. 

Since the mass resolution in the CMS detector depends 
strongly on the pseudorapidity (//) of the reconstructed par- 
ticles, the analysis is performed in two "channels": barrel 
(if both muons have \t]\ < 1.4) and endcap (elsewhere), 
that are then combined for the final result. The determi- 
nation of the signal efficiency in this analysis depends on 
MC simulation, that is validated through two samples of 
fully reconstructed B decays: (1) the decay B ± — > J/if/^, 
that provides a high-statistics sample to allow fine-grained 
comparisons; (2) the decay B s — > J/ift(p tnat is essential to 
compare B s mesons in data and MC simulations and to es- 
timate systematic uncertainties for the analysis efficiency. 

2.1 The signal channel B s(d) -> y^yr 

The events are selected with a two-level trigger system. 
The first hardware level only requires two muon candi- 
dates, while the high-level software trigger uses additional 
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Fig. 2. Comparison of MC signal and sideband data distributions, 
for the transverse momentum (top left), the 3D pointing angle 
(top right), the flight length significance (bottom left), and the 
isolation (bottom right). The MC histograms are normalized to 
the number of events in the data. 
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Fig. 3. Comparison of measured and simulated B± — > J/ i/fX"* dis- 
tributions, for the transverse momentum (top left), the 3D point- 
ing angle (top right), the flight length significance (bottom left), 
and the isolation (bottom right). The MC histograms are normal- 
ized to the number of events in the data. 



2.3 The normalization and control channels 



information from the silicon tracker. For the offline event 
selection, variables related to the muons, the primary ver- 
tex, and the B s candidate with its associated secondary ver- 
tex are calculated. The candidate's secondary vertex and 
its momentum are used to select a matching primary vertex 
based on the distance of closest approach. The flight length 
significance is computed as the ratio of the flight length to 
its error: S3D = hd/o~(hd)- The pointing angle a is defined 
as the angle in three dimensions between the B candidate 
momentum and the vector from the primary vertex to the 
B decay vertex. The isolation / is determined from the B 
candidate transverse momentum and other charged tracks 
in a cone with radius AR - ^J{Arj) 2 + {A<p) 2 = 1 around the 
B momentum as / = PT (B)+£ kPT ■ A comparison of several 
distributions for signal MC and for sideband background 
data events is shown in Fig. [2] 



2.2 Background 

The background sources that mimic the signal topology 
can be grouped into three categories: (1) qq events (where 
q — b, c) with q — > Xpv (prompt or cascade) decays of 
both q hadrons. (2) Events where a true muon is combined 
with a hadron misidentified as a muon (punch-through or 
in-flight decay of a hadron). (3) Rare B°, B ± , Z?° and Ab 
decays, mostly from semileptonic decays: peaking back- 
grounds from rare decays, where a heavy particle decays 
into a pair of hadrons (examples include — > K + K~ , 
Ab — > pK ) and non-peaking backgrounds from rare B°, 
B ± and BP decays, comprising hadronic, semileptonic, and 
radiative decays. 



The reconstruction of B ± — > J/if/K* 1 candidates starts from 
two opposite-charged muons with invariant mass 3.0 < 
m /Jl/j2 < 3.2 GeV, which are combined with a track, as- 
sumed to be a kaon, fulfilling pj > 0.5 GeV. The dis- 
tance of closest approach between all pairs among the three 
tracks is required to be less than 1 mm. All three tracks are 
used in the vertex reconstruction. 

The B s — > J/i//<p decay channel is used to check that 
exclusive B s meson decays are correctly simulated in MC. 
The reconstruction candidates starts from two opposite- 
charged muons, which are combined with two tracks as- 
sumed to be kaons, fulfilling pj > 0.5 GeV and [77I < 2.4. 
The distance of closest approach between all pairs among 
the four tracks is required to be less than 1 mm. All four 
tracks are used in secondary vertex reconstruction. The 
two muons of the candidate must have an invariant mass 
3.0 < < 3.2 GeV. The two kaons must have an in- 

variant mass of 0.995 < m KK < 1.045 GeV and AR < 0.25. 

Comparisons of several distributions for the normaliza- 
tion and control candidates, between MC simulation and 
sideband-subtracted data, are shown in Fig.[3]and Fig. [4] 

2.4 Main systematic uncertainties 

The main source of systematic uncertainty for the branch- 
ing fractions and the production mechanisms is the error 
of fo-quark hadronization fractions f u and f s and it is equal 
to 12.5%. The acceptance dependence on different produc- 
tion mechanisms has been evaluated for signal and nor- 
malization samples and the difference results in 4%. The 
uncertainty on the analysis efficiency is determined as the 
quadratic sum of all the efficiency differences between data 
and MC simulation in the control sample (7.9%). For the 
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Fig. 4. Comparison of measured and simulated B — > 7/(4"0 distri- 
butions, for the transverse momentum (top left), the 3D pointing 
angle (top right), the flight length significance (bottom left), and 
the isolation (bottom right). The MC histograms are normalized 
to the number of events in the data. 

Table 3. Optimized selections for the best upper limit. 



Selection 


Barrel 


Endcap 


Unit 


Pt>, > 


4.5 


4.5 


GeV 


PT.H2 > 


4.0 


4.0 


GeV 


Pt.b > 


6.5 


6.5 


GeV 


I > 


0.75 


0.75 




X/dof < 


1.6 


1.6 




a < 


0.050 


0.025 


rad 


S 3D > 


15.0 


20.0 




4a > 


n/a 


0.015 


cm 



Table 4. Final results for each signal "channel" used for the ex- 
traction of the upper limit. 



Final results 

^(inSB) 

N" p . 

peak 

N obs 

Final results 



B s — > barrel 
0.0036 ± 0.0001 

0.60 ± 0.35 
0.071 ±0.020 
2 



B — > fi + ii barrel 
0.0035 ±0.0001 

0.40 ± 0.23 
0.245 ± 0.056 




B s — > endcap B — > endcap 



A^(inSB) 

peak 

N„ hs 



0.0021 ± 0.0001 

0.80 ± 0.40 
0.044 ±0.011 
1 



0.0019 ±0.0001 

0.53 ± 0.27 
0.158 ±0.039 
1 
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Fig. 5. Unblinded dimuon invariant mass distributions in the bar- 
rel (left) and the endcap (right) channels. 



Including all statistical and systematic uncertainties the 
following upper limits with CLs method are extracted: 



BF(B S 
BF(B 
BF{B S 
BF(B 



yuV) < 1-9 x 1(T 8 at 95% CL 



fi + tT) < 4.6 x 1(T 9 at 95% CL 
(i + tT) < 1.6 x 1(T 8 at 90% CL 
H + H~) < 3.7 x 10~ 9 at 90% CL 



background estimation, candidate counts from the side- 
bands are interpolated in the signal boxes to estimate the 
background candidate yield. This results in a difference of 
4%. For the peaking rare B decays, a set of rare hadronic 
backgrounds was generated, simulated, and passed through 
the full detector simulation and evaluated taking into ac- 
count the CMS muon misidentification. 



3 Results 

To optimize the analysis selection, a grid of selection re- 
quirements was defined and searched for the best expected 
upper limit (Tab. [3}. 

Tab. [4] summarizes all numbers relevant for the extrac- 
tion of the upper limit: the final efficiencies, the number of 
observed background events in the sidebands, the expected 
number of resonant and non-resonant background events 
in the signal windows and the number of observed events. 
Fig.|5]shows the "unblinded" invariant mass distributions. 



4 Conclusions 

A search for the rare decays B s — > yU + /t/~ and Bo — > pfpT 
has been performed on a data sample of pp collisions at 
■\fs = 7 TeV corresponding to an integrated luminosity of 
1.14 fb" 1 collected by CMS. The observed event yields are 
consistent with those expected adding background and SM 
signals. A combination of this result with the LHCb public 
result using 0.34 fb" 1 has been performed 0. The combi- 
nation results in an upper limit on the branching fraction 
of BF(B S -> fi + fT) < 1.1 x 10" 8 at 95% CL. 
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